High temperature superconductivity in dimer array systems 
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Superconductivity in the Hubbard model is studied on a series of lattices in which dimers are 
coupled in various types of arrays. Using fluctuation exchange method and solving the linearized 
Eliashberg equation, the transition temperature T c of these systems is estimated to be much higher 
than that of the Hubbard model on a simple square lattice, which is a model for the high T c cuprates. 
We conclude that these 'dimer array' systems can generally exhibit superconductivity with very high 
T c . Not only d-electron systems, but also p-electron systems may provide various stages for realizing 
the present mechanism. 
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PACS numbers: 74.20-z, 74.20Mn 



Search for high temperature superconductivity is one 
of the most challenging problems in solid state physics. 
Since one of the bottleneck for high T c is the generally 
low energy scale of the phonons, purely electronic pair- 
ing mechanism have attracted much attention from the 
early days. The discovery of high T c superconductivity 
in the cupratestl has certainly boosted interest along this 
line. In particular, the possibility of superconductivity in 
the Hubbard model on a square lattice, a simplest purely 
electronic model for the cuprates, has been intensively 
studied, where spin fluctuation theories using fluctua- 
tion exchange (FLEX) method have estimated T c of the 
d ;c 2_j / 2-wave superconductivity to heJD(0.01t) (t is the 
nearest neighbor hopping integral) .□□ Since t ~ 0.4eV 
for the cuprates, this estimation is consistent with the 
experimentally observed T c of up to 150K, but at the 
same time one should note that it is two orders of mag- 
nitude lower than the kinetic energy scale t. 

In our view, one reason for this reduction of T c is that 
the d x 2_ y 2 gap function changes sign, having nodal lines 
that intersect the Fermi surface. Such a sign change 
in the gap function is a 'necessary evil' for spin fluc- 
tuation pairing in that the pair scattering processes 
[k t, -k 1} -> [k + Q T, —k - Q |], mediated by the 
dominant spin fluctuations with wave vector Q, have 
to accompany a sign change in the gap function <fi, i.e., 
(j>(k)<f>(k + Q) < 0. In the case of the square lattice, 
this requirement leads to the d x 2_ y 2 pairing, where the 
gap function has nodes that intersect the Fermi surface, 
thereby resulting in a reduction of T c . 

From this viewpoint, two of the present authors have 
recently proposed that superconductivity with much 
higher T c can be achieved in systems having two pocket- 
like nested Fermi surfaces, where the dominant pair scat- 
tering processes take place between the two Fermi sur- 
faces, so that the sign change in the gap function takes 
place across, not on, the Fermi surfaces.0 Along this line, 
Hubbard models on a two band lattice (Fig.Ma and a 




FIG. 1. The two band lattice studied in ref. H The thick- 
ness of the lines represent the magnitude of hopping integrals. 



four band ona3 have been studied using FLEX method, 
where high T c has been obtained solving the linearized 
Eliashberg equation. 

If we focus on the two band model shown in Figjl], 
the obtained gap function indeed changes sign across the. 
two bands but stays nearly constant within each band.El 
Such a form of the gap function can be interpreted in 
terms of a real space picture by looking at the system as 
dimers coupled in an array. Namely, it is natural to con- 
sider that the singlet pairs are mainly formed between 
electrons residing on different sites within a dimer due 
to the antiferromagnetic interaction cx t 2 d /U between the 
two sites. Then, the gap function has s wave symme- 
try within each band because the pairing occurs within 
a unit cell, while it changes sign across the two bands 
because the pairs are formed between different sites]] 

This intuitive picture has led us to consider here the 
Hubbard model on a series of lattices in which dimers 
are coupled in various types of arrays. We conclude that 
these 'dimer array' systems can generally exhibit super- 
conductivity with very high T c . 

Let us now start with a dimer array Hubbard system 
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FIG. 2. Dimer array on a square lattice. 



on a square lattice shown in Fig.||.l Here, c\ a / 'a a / 'n c ia and 
d\ a /di„ /nfa- denote the creation/annhilation/number op- 
erators on the lower and upper sites of the i-the dimer, 
respectively. We assume td > t', where td is the hop- 
ping within the dimers, while the hopping t' couples 
the dimers. When [7 = 0, the bonding and antibond- 
ing bands have band width of 8t' each, while the level 
offset between the bands is ltd- So as far as td is not 
too large compared to t', the two bands have a certain 
amount of overlap. Then, if the band filling na is not too 
far away from half- filling (n = 1), both bands cross the 
Fermi level to result in two pocket-like Fermi surfaces, 
one of which lies around k = (0, 0) and the other around 
k = (ir,n) (Fig.||). The two Fermi surfaces are nested 
to some degree when n is close to half-filling, so that 
when U is turned on, antiferromagnetic spin fluctuations 
that mediate inter-Fermi-surface pair scattering (arrows 
in (FigJ|)) arises, which should lead to a superconduc- 
tivity with a gap function that changes sign across the 
Fermi surfaces but not within each Fermi surface. 

On the other hand, td should not be too small, namely, 




FIG. 3. Fermi surfaces of the dimer array system on a 
square lattice. Q denotes the nesting vector, i.e., the peak 
position of the spin susceptibility, while AQ represents the 
spread of the spin susceptibility around Q. 



the overlap of the two bands, and thus the size of the 
Fermi surfaces, should not be too large for high T c su- 
perconductivity to occur. This is because when the Fermi 
surfaces are small compared to the spread AQ of the peak 
in the spin susceptibility, inter-Fermi-surface pair scat- 
terings V[fe T, ~k i] G A -> V[fe + q T, -k - q j] e B 
all have large contribution, while when the Fermi sur- 
faces are large compared to AQ, the pair scattering pro- 
cesses having large contribution will be restricted to a cer- 
tain combination of initial states k 6 A and final states 
k + q e BE 

In order to verify-the above-.expectation, we use the 
multiband versioroEJ of FLEX0, which is kind of a self- 
consistent random phase approximation, to obtain the 
renormalized Green's function of the Hubbard model. 
Then, T c is determined as the temperature at which the 
eigenvalue A of the Eliashberg equation, 
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reaches unity. Here, G is the renormalized Green's func- 
tion matrix (with Z,m,etc... labeling sites in a unit cell) 
obtained by FLEX. V is the spin singlet pairing in- 
teraction matrix given by V(q) — ^V sp (q) — ^V c h(q), 
where the pairing interaction due to spin fluctuations 
(sp) and that due to charge fluctuations (ch) are given 
as V r sp ( c h) = t^ 2 Xsp(ch) using the spin and charge suscep- 
tibilities x Bp (ch)(q) = Xirr(q)[l - (+)E/Xirr(<7)r\ where 
Xirr(<?) is the irreducible susceptibility matrix Xirr(<z) = 
— jj-J2 k G(k + q)G(k) (N is the number of fc-point 
meshes). In the present paper, we take up to 64 x 64 
fc-point meshes and up to 4096 Matsubara frequencies. 

As mentioned above, the FLEX+Eliashberg equation 
approach is known to give reasonable T c of order 0.0 It 
for the single band Hubbard modelcl which is a model 
for the high T c cuprates. The occurrence of supercon- 
ductivity in the Hubbard model contradicts with some of 
the numerical studies, but we believe this to be due to 
finite size effects in the numerical calculations. In fact, 
a quantum Monte Carlo calculation that pays special at- 
tention to the discreteness of the level spacing in finite 
size systems has detected large enhancement of the pair- 
ing correlation function, which ia-at least qualitatively 
consistent with the FLEX results.til 

We have performed FLEX calculation for n = 0.95 
and n = 0.9 with U — 8, varying (t',td) in a certain 
range. The gap function obtained by solving the lin- 
earized Eliashberg equation (Fig.^J) does not change sign 
in each band, while it changes sign across the two bands, 
as expected. As a consequence, T c , shown in Figj^, 
reaches ~ 0.13 for n = 0.95, which is 4 ~ 5 times higher 
than the typical T c of the Hubbard model, jestimated by 
the same way, on a simple square lattice.fjtl 
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FIG. 4. The gap functions for the two bands A and B with 
U = 8, n = 0.95, t' = 0.8, t d = 1.3, and T = 0.15. 




FIG. 5. T c of the dimer array Hubbard model on a square 
lattice as a function of (t',t d ) with U = 8, n = 0.95(a), or 
n = 0.9(b). 



As t' /td increases, the Fermi surfaces become large, re- 
sulting in a reduction of A due to the reason mentioned 
above. At the same time, the antiferromagnetic spin fluc- 
tuations increases rapidly because the nesting between 
the two Fermi surfaces becomes better for large Fermi 
surfaces. Thus, the tendency towards magnetic ordering 
dominates over superconductivity. On the other hand, 
when t' /td is small, only one band intersects the Fermi 
level, so that the Fermi surface nesting is degraded and 
T c becomes low. Consequently, high T c is obtained in a 
certain optimized regime of (t',td) as seen in Fig.||. At 
n = 0.9, T c becomes slightly lower, but still reaches 0.1. 





FIG. 6. Dimer array on triangular (a) and honeycomb (b) 
lattices. 
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FIG. 7. T c of the dimer array systems on triangular (a) or 
honeycomb (b) lattices with U — 8 and n — 0.95. 



We now show that high T c superconductivity in dimer 
array systems is quite general by looking into triangular 
and honeycomb lattices shown in Fig.o(a) and (b), re- 
spectively. T c for n — 0.95 and U = 8 on these lattices 
are shown in Fig.^ as functions of (t',td), which again 
exceed or come close to 0.1. These results suggest that 
dimer array systems can generally exhibit high T c super- 
conductivity on various types of lattices. 




FIG. 8. (a)Dimer array system (consisting of p w orbitals 
here) on a honeycomb lattice in which the site energy of the 
unhatched orbitals is higher by e than that of the hatched 
ones. (b)The effective lattice of (a). 



example, let us consider a dimer array system on a hon- 
eycomb lattice where the site,energy differs by e between 
neighboring sites (Fig.|](a)) .Ej At quarter filling (n = 0.5) 
and in the limit of large e, the system becomes equiva- 
lent to a half-filled dimer array system on a triangular 
lattice with an effective hopping t' cS = t 2 je between the 
dimers (Fig.||(b)). Consequently, we can expect high T c 
when the value of s is tuned so that U /t' cS is optimized. 
In Fig.|, we show T c for U = 2, t' = 1, t d = 0.3B and 
n = 0.5 as a function of e. T c is found to exceed 0.02, 
which again corresponds to few hundred K if t' = 1 cor- 
responds to few eV as in 2p systems. 

To summarize, we have shown that 'dimer array' Hub- 
bard systems can generally exhibit superconductivity 
with very high T c . Not only d-electron systems, but also 
p-electron systems may provide various stages for realiz- 
ing the present mechanism. 

We would like to acknowledge H. Aoki, J. Akimitsu, 
and Y. Zenitani for valuable discussions. Numerical cal- 
culations where performed at the Supercomputer Center, 
Institute for Solid Stated Physics, University of Tokyo. 
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FIG. 9. T c of the system shown in Fig|8j(a) as a function of 
e with U = 2, t' = 1, t d = 0.5 and n = 0.5. 



So far we have focused on parameter regimes with large 
U/t' and U ltd, where we have d-electront3 systems in 
mind. Namely, if t' and td are few hundred meV, which 
is typical for ci-electron systems, T c ~ O.lt' corresponds 
to a very high temperature of few hundred K. Then our 
next question is : can we obtain similarly high T c in sys- 
tems with smaller U/t', as in p-electron systems ? In 
2p-electron systems such as B, C, and/or N compounds, 
the on site Coulomb repulsion and the hopping integrals 
are both typically few eV, so that U/t ~ 0(1). 

In such a case, we have found that T c estimated by 
FLEX+Eliashberg equation is very low (or does not ex- 
ist) not only on conventional lattices, but also on the 
lattices considered above as well. Nevertheless, we now 
show that high T c can be achieved even for systems with 
small U/t' by introducing a level offset between neighbor- 
ing sites, thus reducing the effective band width. As an 
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